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Abstract 

Complex structures are determined for surfaces with S 2 and T 2 topologies generated by the 
dynamical triangulation method. For a surface with 5* topology the spacial distribution of the 
conformal mode is obtained, while for the case of T 2 topology the distribution of the moduli parameter 
is calculated. It is also shown that the network of Feynman diagrams of massive <j> scalar theory 
has a unique complex structure. This gives a numerical justification of the hadronic string model for 
explaining the n-particle dual amplitude. 



'Based on a contribution paper for the proceedings of the international conference "Frontiers in Quantum Field Theory" 
celebrating Professor Kikkawa's 60th birthday, held at Osaka University, Osaka, Japan, December 14-17, 1995. 
Presented by T.Yukawa. 

^E-mail address: kawaih@theory.kek.jp 

*A JSPS Research Fellow, E-mail address: ntsuda@theory.kek.jp 
§ E-mail address: yukawa@theory.kek.jp 



2.0 3.0 4.0 5.0 6.0 2.0 3.0 4.0 5.0 6.0 2.0 3.0 4.0 5.0 6.0 
r r r 



Figure 1: Resistivities for the case of pure gravity. The size of (a), (b) and (c) are 4000, 8000 and 16000 
triangles respectively. 

1 Complex structure in the Polyakov string 

For quantizing the Polyakov string one usually considers a set of closed string world sheets of the 2- 
dimensional Euclidean space-time [[l] . Each world sheet is considered to be a closed and orientable manifold 
with various topology. Once the Riemannian metric g a b is given on this manifold, it is always possible to 
introduce a complex structure on the surface. For example, the Riemannian metric of a surface with S 2 
topology can be written by a suitable choice of local coordinates as 

.9 = g a bdx a dx h = e a( -^dzdz, 

where z = x + iy, and e a ^ is the conformal factor. Therefore the problem of quantizing two-dimensional 
gravity is reduced to considering the quantum fluctuations of the conformal mode and the complex 
structure. This is the basic assumption one usually employs in the Liouvillc field theory of non-critical 
Polyakov string. 

Now, we would like to examine this assumption constructively, by determining the complex structure 
for those surfaces obtained numerically by the dynamical triangulation(DT). The Monte Carlo simulation 
of the random surface by the DT method Q has been regarded as a numerical exercise of the matrix 
model ]3| which is known to give the correct critical exponent identical to the continuous field theory, i.e. 
the Liouville field theory Q]. In what follows we shall show that the DT method also exhibits its explicit 
correspondence to the continuous field theory in the sense that those surfaces created by DT numerically 
can be decomposed to a unique complex structure and conformal modeQ. 

2 The basic idea -resistivity- 

Our method is based on the observation that the resistance of a homogeneous conducting sheet is in- 
variant under the local scale transformation. This can be shown by considering the resistance of a small 
rectangular section with length a and width b of a conducting sheet. The resistance between two sides of 
length a is given by 

„ a, 
R = r-, 

where r is the resistivity constant. The resistance R is apparently invariant under the scale change 
a t— » arj, b i— > bij. Therefore, the local scale factor does not affect the resistance. 

Experimentally, the resistance measurement of a conducting sheet is carried out by picking up four 
points on the surface which we specify by four complex numbers {zi}. With a source of current / placed 
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Figure 2: Resistivities for the case of 2 scalars. 



at z 3 and a sink of the current at 2:4 the voltage drop from a point z\ to a point z 2 is written as 

V(z 1 ) - V(z 2 ) = R(zi,z 2 ; z 3 , z 4 )I. 

In the case of surface with S 2 topology it can be regarded as an infinite flat sheet and we obtain from 
Gauss' law 

T 

R = -— \n\[z 1 ,z 2 ;z 3 ,z i ]\, 

where 

r 1 Zl — Z3 Z 2 — Z4 

[zi, z 2 ; Z3, Zi\ — 

z 2 — 23 z \ — z i 

is known as the anharmonic ratio. 

We regard the dual graph of a triangulated surface as a trivalent network of registers each having lfL 
We pick up two vertices for applying current with 1A, and solve the Kirchhoff equation by the Jacobi 
iteration method. Then, we measure voltage drops between other pair of vertices. Using the SL(2, C) 
transformation, 

az + b, » 
z 1 > -(ad — bc = I), 

cz + a 

invariancc property of the anharmonic ratio we can fix three points among {zi\ at z 3 = 0,24 = 00 and 
z 2 = \ without changing the resistance by an appropriate choice of four complex parameters (a, b, c, d). In 
practice, among 4! combinations of measuring arrangements there are known to be only two independent 
measurements, while we have three unknowns, namely one complex z\ and one real r. In order to fix 
them unambiguously we add one more measuring point z 5 . 



3 Complex structure of S* 2 surfaces 

Let us show results of the resistivity measurement for random surfaces with S 2 topology simulated by the 
DT method with N triangles. If a conducting surface is uniform and homogeneous, we should get a fixed 
resistivity regardless of points of measurements. Fig. I shows the distributions of resistivity constant for 
surfaces of three area sizes with no matter fields coupled(the pure gravity). The distributions peak at 
about r = 2.6 and they get sharper as the size grows. Broader distributions in small surface simulations 
reflect the mesh structures of circuit networks, and we can expect that the peak will eventually become 
the 5-function type distribution as the area get infinity. On the contrary, when surfaces are coupled with 
matter fields with the central charge c bigger than 1 , peaks of the resistivity distribution tend broader as 
the size grows. 
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Figure 3: Resistivities for the case of 3 scalars. 
3 Ising spins 
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Figure 4: We plot resistivities for 3 Ising spins both case of off-critical (a) and near critical (b). 



Fig.2(c = 2) and Fig.3(c = 3) show the tendency clearly. This phenomenon suggests the absence of 
a unique continuous limit for cases c > 1 as expected from analytic theories. The c = 1 barrier can be 
seen more significantly in the measurement of a surface with three Ising spins coupled(Fig.4). It should 
behave as a system with the c = | matter coupled at the critical point, while it is expected to be in the 
same universality class as the pure gravity off the critical point. Indeed, our measurement shows this 
behavior. Broadening of the resistivity for bigger size naturally implies the transition of the surface to 
the branched polymer phase beyond c = 1 . 



4 Distribution of the conformal factors 

Once we obtain the resistivity r for a surface, we can locate the position of each triangle by two mea- 
surements of resistances, for example 

T 

i?(zi,l;0,oo) = -—In N, 

and 

i?(z l7 0;l,oo) = --L]n|l-*i|. 

By sweeping all the rest of triangles with fixing three points which correspond to (1,0, oo), we can 
determine the position of all the triangles of this surface. The point density obviously represents the 
distribution of conformal factors proportional to 

lim |z 4 | 4 < e «*(0) e aa(l) e a*(z 4 ) e a<T(*i) > _ 

24^00 
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Figure 5: A plot of the obtained coordinates of all vertices for c — with five independent configurations 
superimposed. As a matter of convenience all points are mapped into a circle with the center at (51 2 ) 
and radius 



In Fig. 5 we show the distribution of triangles whose position z is mapped within a circle by the transfor- 
mation 



z + a — 1 



with a = §+*2 • There exists a theoretical prediction of the distribution^, but it has not been presented 
in a form ready for the direct comparison to the numerical experiments. 



5 Complex structure of the T 2 surface 

As far as two dimensional orientable manifolds of S 2 topology are concerned they all have the same 
complex structure, i.e. any manifold can be transformed to the other by a combination of a general 
coordinate transformation(Di/Jeo) and the Weyl transformation ( Weyl). On the other hand for the surface 
with T 2 topology the moduli space 

{gab} /Weyl <g> Diffeo 

is spanned by a complex plane of moduli t. It is defined by the ratio of the integrals of an Abelian 
differential 

where integration contours (we call them as the a-cycle and the b-cycle) are chosen to be two closed 
paths on the manifold which cross each other only once. Here, j^dx^ is a harmonic 1-form and j^dx^ 
is its dual. If we regard as a current density on the torus, r § a j^dx >l and § a j^dx^ correspond to the 
voltage drop along the a-cycle and the total current crossing the a-cycle, respectively. 

In practice, we select two closed paths made up by connecting edges of triangles, which intersect only 
once. Then we cut the surface along one of the path for which we choose the b-cycle, and apply constant 
voltages(lV) in between neighbouring triangles dual to the cut(Fig.6). Solving the Kirchhoff equation we 
can determine all the currents and j M along links between neighbouring triangles in the dual graph. 
By construction of the network the following two integrals are trivial; 

j> j^dx^ = ji judxi* = 0. 
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Figure 6: 



Thus r is obtained by measuring total currents crossing two cycles; 

_ ijb 

where I a and h represent the total currents crossing the a-cycle and the b-cycle, respectively. Instead of 
measuring the resistivity of this surface we employ the value determined previously for the surface with 
S 2 topology and corresponding central charge. 

Fig. 7 shows the distribution of r for configurations of 4,000 and 8,000 triangles. Each value of t is 
transformed appropriately by the SL(2, Z) transformation 

t i ► ^— (ad — be = 1) 

cr + d y ' 

so that it is in the fundamental domain. By summing over real t's with a fixed imaginary part we get the 
distribution for surfaces with two, four and eight thousand triangles in Fig. 8, which should be compared 
with the theoretical prediction 
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n=l 



The Liouvillc theory predicts for a surface with c scalar fields 



-2(c-l) 



i.e. for simulations with c > 1 matter the surface will be unstable to the polymer like torus configuration 
in the limit Imr — > oo. 



6 Hadronic strings -a root of the string model- 

Thc n-particle amplitudes of the massive </> 3 scalar theory corresponding to the diagram G shown by 
Fig. 9 is given by 

i G ( P1 ,p 2 ,..., Pn )= / nAin^nw 

i=l i * a=l 
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Figure 7: Plot of the moduli(r) on the complex-plane mapped to the fundamental domain with a total 
number of triangles of 4, 000 and 8, 000. Each dot corresponds to the modulir for a sample of DT-surface. 
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Figure 8: Four density-distributions (histogram) of imaginary part of moduli Ti for the case of pure- 
gravity with the topology of the torus. These are obtained by integrating over of the real part of moduli 
ti on the complex plane. A scale of vertical line is arbitrary. 
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Figure 9: 

where P a is the net momentum flowing into the vertex a. Here, Nq stands for the total number of vertices, 
and N% for the total number of internal lines. By the Laplace transform of a propagator, 

1 r°° 

qf + m 2 J a 

the amplitude is expressed after integrated over {qi} as 

i a i 

where a(i) and b(i) are two vertices connected by the propagator i. This equation has been often noticed 
on its resemblance to the electric circuit 0, i.e. 

Ui — » resister,a; a — ► voltage, p a — > current. 

Then, Q = g»(»^f6W — jg regarded as the total heat generated by the circuit. It has been conjectured 
that as the mesh becomes finer the surface will tend to be a uniform, homogeneous and continuous 
conductor. In this case after averaging over all n-particle diagrams we will obtain 

< I >~ e "^»." flab(a)PaPb , 
where a is an appropriate resistivity of the surface and the resistance is same as that of the flat conductor; 

Rab oc In \z a - z b \ 

We have made the measurement of resistivity for networks with S 2 topology constructed by the same 
DT method as described previously. This time each link has resistance distributed randomly with the 
probability 

Here, we have chosen the parameters to be D = and m = 1 in the simulation. Fig. 10 shows the size 
dependence of the resistivity distribution. As the size gets larger the distribution becomes sharper similar 
to the case of the pure gravity. This gives a numerical support of the theoretical conjecture. 
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Figure 10: Resistivities for the case of pure gravity when each resistance is randomly distributed. 
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